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1. Introduction 

Over the last year and a half, revolutionary progress in understanding the dynamics 
of four-dimensional supersymmetric gauge theories has been made by Seiberg and collab- 
orators . One spectacular result is the exact low-energy effective Lagrangian for Af = 2 
supersymmetric SU (2) gauge theory obtained by Seiberg and Witten || . The Af = 2 vec- 
tor multiplet contains an adjoint scalar whose non-zero vacuum expectation value (vev) 
breaks SU(2) to U(l). For large vev compared to the scale A set by the gauge coupling, 
one can write an effective Lagrangian in terms of a U(l) gauge multiplet. For small vev, 
in the past we would have expected SU(2) color confinement and a very different (and 
inaccessible) description. 

Surprisingly, it turned out that the effective theory is a U(l) gauge theory for any 
expectation value, and the naive unbroken SU{2) regime is not present at all. Instead, 
two singularities of the effective action appear in the strong coupling regime. At one of 
these, the monopole visible in the semiclassical treatment becomes arbitrarily light. The 
effective Lagrangian is again a U(l) gauge theory but now written in terms of a vector 
multiplet containing the dual (magnetic) gauge field with the standard local coupling to 
the monopole, as well as a hypermultiplet describing the monopoles. The other singularity 
is isomorphic to this, with the role of the monopole taken by a charge (1,1) dyon. 

Thus the Af = 2 theory does not confine. Surprising as this result may be, it does not 
drastically contradict previous expectations, mostly because of the presence of the massless 
scalar. Theories which do not have such a scalar and are widely believed to be similar to 
pure (bosonic) gauge theory are Af = 1 SYM gauge theory as well as supersymmetric 
QCD. These theories should confine, and Seiberg and Witten showed that this could be 
explained in Af = 1 theories obtained by adding an Af = 2 breaking mass term: it is the 
result of monopole condensation. 

These results have been extended by Argyres and Faraggi [|] and by Klemm, Lerche, 
Theisen and Yankielowicz |4]] to pure SU(N) gauge theory.* The Af = 2 theory now has 
N — 1 moduli (say, the eigenvalues of the adjoint scalar vev (0)) and in the semiclassical 
regime the gauge symmetry is broken to [/(l)^ -1 . Each factor contains a monopole 
solution and by varying the moduli one can drive any of them massless. There are A" 
points in moduli space at which N — 1 monopoles becomes massless, and these points 
become the ground states of the associated Af — 1 theory. 

* It should be noted that the precise dependence of the effective action on the moduli proposed 
in equations (4) and (14) of the second paper in j|] is ambiguous for general N. We follow the 
unambiguous results in equations (6) and (9) of || . 
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In this paper we continue the discussion of the physics of these theories begun in 
0,|3],|]]. Perhaps the most striking result is the following. In the weakly coupled regime, 
the theory is a ^(l)^ -1 gauge theory with a discrete gauge symmetry Sn permuting the 
£7(1) factors. This discrete gauge symmetry is spontaneously broken by the Higgs vev - 
for example, at a generic point in moduli space, every charged multiplet has a distinct 
mass. It turns out that this is true everywhere in moduli space, even at the vacua of the 
TV = 1 theory. This leads to a non-trivial spectrum of light massive particles in the M = 1 
theory and rather surprisingly, a spectrum of distinct string tensions in the different U(l) 
factors. Since the theory contains particles with charges in any pair of U(l) factors, the 
limiting (infinite distance) string tension is the lowest of these, but mesons and baryons 
bound with the higher string tensions will exist as sharp resonances for sufficiently small 
M = 2 breaking. 

We also consider the large iV limit of the theory and compare with expectations from 
previous work. Let us first say a few words about possible large N limits. The bare coupling 
constant is always taken to zero as g$ ~ 1/N to get a theory with a planar diagram 
expansion. Now we will discuss only the leading (two derivative) effective Lagrangian, 
and in M = 2 SYM this receives no perturbative contributions beyond one loop. Since 
instantons are suppressed as e~ N ^ 9 ° , it may at first sound like there is little to do. However, 
this statement turns out to be naive. 

In the M = 2 theory, there are two natural regimes to consider. To give all charged 
bosons masses of at least O(7V ), the difference between any pair of eigenvalues of <fi must 
be at least O(iV ), and the typical difference will be O(TV).* The density of states will be 
O(N). We refer to this as the naive semiclassical regime as monopole masses are 0(7V), 
and instanton corrections do not survive in this limit. The cause of this is not only the 
e ~ N /9 suppression. Rather it is a non-trivial consequence of the need to introduce iV 
distinct mass scales. 

To make instanton corrections survive, one must take the spacing between eigenvalues 
to be 0(1/N). In this scaling regime, monopoles are light. The vacua of the M = 1 theory 
are of this form, and one can smoothly interpolate between them and a limit in which the 
spacings are A/TV with A becoming large, in which the semiclassical treatment is valid. In 
this sense, there is no large iV transition in the theory. 

Since the effective theory is a spontaneously broken abelian gauge theory, the standard 
expectations of the large N limit - a finite mass gap, O(l) degeneracies in the particle 

* We assume for simplicity that all eigenvalues are real. Otherwise there are configurations 
where the typical difference will be 0{\fN). 
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spectrum, and an effective (p 3 coupling of order 1/N - are not at all obvious. We find that 
at least the last of these is violated near the massless monopole point. 

The traditional definition of the large N limit takes N — > oo before any other limits 
and in particular before the infinite volume limit. Our low energy effective Lagrangian 
should be accurate for processes involving momenta \p\ << A e jy, where the scale is set 
both by the masses of particles we have integrated out and by the couplings of higher 
derivative terms we have dropped. Electrically charged particles must be integrated out, 
and it will turn out that the lightest of these has mass ~ A/iV 2 , which will make the 
interpretation of the limit subtle. 



(2.1) 



2. SU(N) supersymmetric gauge theory 

The Af = 2-supersymmetric bare Lagrangian is (in Af = 1 superfield notation) 

ro = Ani/g 2 + 6/2n is the bare gauge coupling. There are two reasons for the explicit 
N dependence. Of course, it is the appropriate one to weigh Feynman diagrams of Euler 
character \ as N x : and thus we will take the large N limit with To fixed. It is convenient 
even at finite N, since it cancels the explicit N in the one-loop beta function, and thus the 
dynamical scale A ~ exp — Imro (at which the running coupling attains a prescribed 0(1) 
value) will have no N dependence. 

The Af = 2 theory has an N — 1 complex dimensional moduli space A4 of vacua. 
In the classical theory these are parameterized by the invariant expectation values Tr <p n 
constructed from the scalar component of A. <fi must satisfy the D-flatness condition 
[<p, <p + ] = and thus we can diagonalize it, and use as coordinates for A4 the eigenvalues 
4> m with permutations identified. At a generic point <fi breaks the gauge symmetry to 
t/^l)^ -1 and a low-energy effective Lagrangian can be written in terms of multiplets 
(Ai, Wi). We will use a l U{l) N ' notation in which 1 < i < N and £\ A t = 0. We denote 
the scalar component of Ai by a^. 

The Af = 2 effective Lagrangian is determined by an analytic prepotential T and takes 
the form 

(2.2) 



£e// = Im 4^ 



J d 4 9 d l T{A)A l + 1 J d 2 9 d i d j J r (A)W i W j 



In the classical theory, JF cZ (A) = 4p — Aj/N) 2 . For large <j) m the gauge coupling 

is weak at the scale of symmetry breaking and a good approximation to T would be 
obtained by adding the one-loop quantum correction: 

fi^B^-VK^ (2-3) 
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One can renormalize and define A to absorb T c i into this expression.* In fact, M = 
2 supersymmetry forbids further perturbative corrections |J . This is not to say that 
perturbation theory is trivial but that it will only produce higher derivative terms in the 
effective Lagrangian. 

The reduced (or BPS saturated) multiplets and their mass formula play a central 
role in the story. The [/(l)^ -1 theory has a lattice of allowed electric and magnetic 
charges, which we write qi and hi, again with £\ q$ = hj, = 0. The charges of vector 
bosons are the vectors q v = (0, . . . , 0, +1, 0, . . . , 0, —1, 0, . . .). The fundamental representa- 
tion (we use the convenient name 'quark') are q q = (0, . . . , 0, 1, . . . 0) — (1/N, . . . , 1/N). 
The theory contains 't Hooft-Polyakov monopoles as well, whose magnetic charges must 
satisfy the DSZ condition q^ ■ - qW ■ e Z. If we order the eigenvalues 
<pi > 4>i+i, we expect the stable monopoles to be those with charges in successive fac- 
tors: h = (0, . . . , 0, +1, —1, 0, . . .), so we introduce the basis vectors h l m = 5 l m — 8 l m+x . The 
mass of a reduced multiplet is determined by its M = 2 central charge, which is determined 
by its electric and magnetic charges to be 

M = V2\Z\ = V2\a-q + a D -h\ (2.4) 

with cioi = dF/dcii, a result motivated by duality as we discuss below. 

The an associated with individual monopoles are a Dm = h l m aDi, whose inverse is 
&Di = dD,i=i — J2 m <i a Dm- We will always use the indices ij versus mn to distinguish the 
two bases. Finally, we introduce a n = J2i< n a i w hh inverse = q™a n = a n= i — a n= i-i. 
These are canonically conjugate to the aDm in the sense that a n = —OJ-'d/ 'dciDn- 

We first discuss the naive semiclassical regime and its large iV limit. To break SU(N) 
to t/^l)^ -1 and give every charged multiplet a mass of order iV° or greater, we need to 
choose so that the minimum difference between eigenvalues is O(N ). A representative 
choice is (pij = v(i — (N + l)/2)5ij, where v is the characteristic scale of the vev. The 
classical prediction for the mass of the multiplet Aij is then v\i — j\ and we have a linearly 
rising spectrum with multiplicity O(N) at each mass level. Monopoles are expected to 
have masses AwNy/g^ and do, although it may be interesting to note that the formula 
dDi = d{T c i + J~i)/da,i produces unusual corrections to this, e.g. with log N dependence. 

We now turn to the exact solution . ^(A) is determined in the quantum 

theory by combining analyticity with a physical ansatz for the number and types of its 
singularities, at points where BPS saturated states become massless. What makes it 

* We introduced the extra 'e 3 ' compared with [Q] to be consistent with an extra ^ in the curve 
(2.6) below. It will turn out that this simplifies the formulas at the massless monopole point. 
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possible to study the strong coupling regime and combine the information from different 
regimes is the existence of exact duality transformations on the effective Lagrangian, and 
the generalization of the Witten effect: encircling any singularity in moduli space produces 
a non-trivial Sp(2N—2; Z) transformation on the electric and magnetic charges of all states. 

J- {A) is not a single- valued function of A (as is already clear from ( |2.3|) ) and thus we 
must distinguish the coordinates on M. (for which we retain the name (pi) from the scalar 
components (a^,^) of the fields in the effective Lagrangian. It turns out that J-'(A) is 
most simply expressed in terms of an auxiliary Riemann surface C which varies on M., 
defined by the curve 

.2 _ A 27V 



y z = P(xY - A' 

1 ^ /..,>_ :I TT,.. ... , (2.5) 

2 



P(x) = idet (x -((/))) = ll[(x 



(We generally set A = 1 in the following). The (aoi,(ij) are then integrals of the mero- 
morphic form A = (l/2ivi)(x/y)dP(x) over a basis of one-cycles with the intersection form 
hi ■ qj. 

To define the one-cycles, order the branch points Xi. Let 7$ for 1 < i < N encircle the 
branch points X2i-i and X21, and a m for 1 < m < N encircle the cut running from X2m to 
X2m+i- The 7i are not independent but satisfy Yliili = 0' while the a m are independent. 
Their intersection matrix is (a m , jj) = 8 m j — 8 m j + i and thus we can associate the quarks 
with the cycles 7$ and the monopoles with a m . We also define a set of cycles conjugate to 
the 

We then write 



a-Dm = f A a n = (h A. (2.7) 

The strong coupling regime is controlled by points where monopoles and dyons become 
massless: h ■ do + q ■ a = 0. The vacua of the Af = 1 theory will come from points in 
moduli space at which monopoles coupling to each U(l) become massless. There are iV 
such points, with a simultaneous degeneration of all the a-cycles of the quantum curve. In 
the realization y 2 = P(x) 2 — 1 this will happen when the N cuts are lined up, each with 
one branch point coinciding with the next. In other words, we require P(x) 2 — 1 to have 
N — 1 double zeros and two single zeros. This condition can be satisfied using Chebyshev 
polynomials: 

P{x) = — det (x — ((f))) = T/v ^— J = cos ^iVarccos • 



Pi.rV I = [ — -1 ] U N ^ 



x 2 \ ( x^ (2i 
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X 




We obtain N — l more solutions by complex rotations x — > e l7rr / N x. Since each is associated 
with a ground state of the M = 1 theory, whose Witten index is N, we do not expect other 
solutions to exist.* 



* The full story is subtle as there are partial degenerations of the curve at which N—l periods 
of A vanish. These are limits in which more than two branch points coalesce, and an example 
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The eigenvalues of <j> are non-degenerate, <p n = 2cos7r(n — \)/N. They have spacing 
0(1/N) in the center of the "band" and 0(1/N 2 ) at the edges. The double branch points 
of the curve are <j) n = 2 cosim/N for 1 < n < N and the single branch points are at ±2. 

Let us compute the periods of the maximally degenerate curve Co- Since the cuts 
are all on the real axis, the a periods will be imaginary while the 7 periods will be real. 
Changing variables from x = 2 cos 9 to 9, we have 



P(x) 

y 

A 



cos N9 
% sin N9 
1 x dP(x) 



d(j>t 



2ni y dx 
N 

— cos OdO 

TV 

1 ldP(x 



dx 



(2.9) 



2ni y d<j>i 
1 1 



dx +d(. ..) 



7T X — 0, 



cot N6 sin 9 d6. 



The aD m are integrals around the a cycles. These degenerate to integrals around 
the 9 m . Since A is non-singular, all a Dm = for C - When we vary the curve, we could 
get two types of contributions. First, the zeroes of P 2 — 1 split. This will be important 
below but here we simply inflate the contour to enclose the new branch points. Second, 
the derivatives of A have poles: 



B r 



.da 



Dm 



34>i 



de 



1 



cot N6 sin 9 



x - (pi 
cos N9 sin 9 



sin N9 



4>i - 4>r. 

1 sin 9 m 

N cos Oj — cos 9 r 



8—0 ri 



(2.10) 



is P(x) = 7}X N — A N . A coalescence of k branch points will cause the periods of A on cycles 
surrounding any two branch points to vanish simultaneously. One can choose a basis of k — 1 
such cycles, and they have a non-zero intersection form, meaning that the associated massless 
particles would not be mutually local. One can show that at such points at least one U(l) will 
not couple to any of the massless particles and so will not confine when an Af = 2-breaking mass 
perturbation forces the massless particles to condense. Thus these are not candidates for J\f = 1 
vacua, where full confinement is expected. We thank P. Argyres, W. Lerche and L. Randall for 
discussions on this point. 
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The matrix B mi is simple in the following basis (appendix A): 



E 1 Tvkm nk(i — i) 
BrmSin-^- = cos — . (Al) 



m=l 



The a m are integrals around the (3 cycles. We find 

/dm 
A 

^ Jo 

2N nm 
= sin — — . 

7T N 

These determine the masses of BPS saturated electrically charged states. Quark 
masses would be given by the 7 periods, which are differences of (3 periods: rrij = 
2v ^ JV (sin ^ — sin 7r ^ 1 ^ ). For odd A, m^+i)/2 vanishes, but since there are no quarks in 
the theory at hand, this is not a difficulty. The 'W bosons' which are present have masses 
rriij = mi — rrij, which are all non-zero at finite A. 

The first comment to make is that these masses are different for particles with charges 
in different U(l) factors. As we commented in the introduction, this is generically true 
at weak coupling: an explicit choice of Higgs vev will break the discrete gauge symmetry 
relating the factors. What may be surprising is that this persists for small Higgs vev and 
at the vacuum relevant for the related J\f = 1 gauge theory. There is simply no symmetry 
under permuting the periods of the curve Co- 

Our second comment is that the mass of the lightest W boson goes to zero in the 
large A limit as mi2 = Atc 2 /N 2 . This is important because it determines the energy scale 
at which our effective Lagrangian breaks down, and we discuss this point below. 

The derivatives of the 

da m 1 f sin6 

Ami = -TT- = 7T ? d9 a T COt N9 

Od>i 2-7T la COS — COS Vi 

(2.12) 

1 r m/N in sine 
= - dO cot N9 

7T Jq COS — COS &i 

This is log divergent, as it should be on physical grounds. The dual gauge coupling in 
each U(l) factor has a positive (IR free) beta function produced by a one-loop diagram 
involving its light monopole. Away from the massless monopole point, it will be cut 



8 



off at the mass of the monopole to produce the same result in each factor as in ||:* 
4:7V /e 2 D ~ r^ n ~ -(i/2Tv)d mn loga Drn . 

To get a finite result one must perturb the curve slightly to give the monopoles small 
masses. This again splits the double zeroes and modifies A. For calculating the coefficient of 
the logarithm, the modification of A does not matter, and the result is given by integrating 
(d\/d(J)i)c Q up to the branch point, which we parametrize as = m — e m /N. The integral 
then simply gives an endpoint divergence, 

e=irm/N-e m /N 



1 sin# , AT/1 

Ami 77 7 7 log Sin NO 

7TN COS Oi — COS 

~ loge m B mi . 

7T 



?=o (2.13) 



The e m are computed by varying the equation y 2 = P(x) 2 — 1. Since we are splitting 
double zeroes we will find e ~ (Sep) 1 / 2 . Using (|2.9j ) we have 

= (^) 2 + 2P^— Sfr 

BP 

= - sin 2 N50 - 2 cos 2 NO V -^Sfa 

-2^^ 5<j>i (2.14) 



e 2 

m 



i 

N 
sin m 

N 
sin m 



i 



using ( |2.10| ). 

Thus the period matrix diverges as 

da. 



D 

mn da Dr 



i 



— d mn log— 

ZTC Am. 



(2.15) 



with A m = Asin^m/A^.* This checks with the expectations for the beta function, and 
confirms the existence of one monopole in each factor (there is no overall N) . 



* The i/ir of ^ is in different charge conventions. 

* The preceding calculation of A m i left out constant terms coming from the bulk of the 
integral, which could change A m . In section 5 we will do this more carefully and show that 
important constants are present, but are compatible with the definition of A m we give here. 
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At this order the different U(l) factors are completely decoupled. Note that the 
decoupled factors are different at weak and at strong coupling (r is diagonal in a different 
basis). Physically, the two limits are controlled by different light degrees of freedom. The 
m'th monopole beta function turns on at the scale A m , and at sufficiently low energy, their 
contribution 5 mn /e 2 D will dominate any structure in r from higher energies, and decouple 
the factors. 

In general, an extended object of size L does not contribute to the beta function 
at energies E > 1/L, so we can interpret the scale A m as an indication of the size of the 
monopole. Semiclassically the monopole size would have been L ~ l/rriw and even though 
here gauge couplings are O(l), we still have A m = m m)m+ i/2v^27r with m mjm +i to good 
accuracy the lightest W mass coupling to that factor. 

The effective Lagrangian around the maximal degeneration is 



The U(l) factors will couple at higher orders in an expansion of the kinetic term, but we 
expect these interactions to be suppressed by powers of p/A m . One physical source of 
this coupling is loops of massive W particles charged under more than one U(l), and such 
interactions will be suppressed by the mass of the Ws. 

An interesting question for contact with the standard large N limit (to be discussed 
below) is whether the interactions are suppressed in the standard way, as <p 2+k /N k . It is 
obviously not so for the superpotential terms in ( |2.16|) , which are completely determined 
by J\f = 2. Nor is there an obvious reason for it to be true of the couplings between U(l) 
factors. 

A M = 1 theory can be obtained by adding the perturbation W = NmTr A 2 to 
the bare superpotential. Following 0] we interpret this in the effective Lagrangian as the 
perturbation W = NmJ2 i (f) 2 , the observable (single- valued on moduli space) equivalent 
to A^mTr A 2 in the semiclassical limit, and thus the exact superpotential is 




(2.16) 




rn 




(2.17) 



rn 
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To calculate its effect we need to change variables from fa to a Dm- We find that the 
vacua W = satisfy 

V2(MM) n = 2Nm^(B- 1 )m(f>i 

(2.18) 

J2 B m{MM) n = V2Nmfa 

n 

and ciDn = 0. This means of course that the kinetic term cannot literally be given by 
( |2.15| ) . Physically, the monopole loop integrals are now cut off by masses produced by the 

~ 1/2 

M = 1 part of the superpotential. We use the prescription of taking a,Dn = (MM) n in 
the M — 2 kinetic term to account for this. 
The equation ( |2.18j ) has the solution 

(MM) n = 2 V2Nm sin ^ (2.19) 

and all of the U(l) gauge symmetries are spontaneously broken. We did not find an 
explicit discussion of the particular M = 2 abelian Higgs theory which appears here in the 
literature, but it is easy to work out* that the particles in a given U(l) factor all have the 
same mass 

m 2 n = 2e 2 Dn (MM) n (2.20) 
and the string tension in the factor is 

n n = 2ir(MM) n . (2.21) 

It is different in each factor, and roughly proportional to the lightest W mass in the factor, 
K n oc mN 2 m n m+i- The gap is non- vanishing in the large iV limit. Because the scalar and 
vector retain the same mass after M = 2 breaking, there is no long range potential between 
two strings (they are 'neutrally stable'). 

Using our prescription to determine the kinetic term, 

An 1 mA nn. 

iir-s 108 ^ 8111 ^ 



1 , . mN 
i- lo g( 



3 (2-22) 



47T A sin ^ 

and there is a weak (logarithmic) variation between the couplings in the different U(l) 
factors. 



* One can simplify one's life by observing that the linear part of the J\f = 2 breaking term in 
W is a pure F term, which using J\f = 2 can be rotated into a D term. This turns the vev into 
(M 2 ) = 2v / 2iVmsin ^ and (M) = and the string solution becomes that of the usual J\f = 1 
abelian Higgs model. 
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3. Physics at finite N. 

The effective Lagrangian ( [2.16|) with the superpotential ( |2.17|) provides an expiicit re- 
aiization of the abelian monopole condensation model of confinement in an SU(N) gauge 
theory. Such models of confinement were much discussed previously, as in 0, on a qual- 
itative level. However, what in retrospect seems an obvious question was not discussed 
to our knowledge: namely, how can an analysis in which the theory looks like a broken 
U(l^N-i g au g e theory avoid having N — 1 distinct types of flux tubes and a corresponding 
multiplicity in the spectrum? 

In fact we see from (|2.21|) that the U(l) factors have differing string tensions. (There 
is a Z2 symmetry a« — > — a/v+i— i, (MM) n — > (MM) N-n which reduces the number of 
distinct string tensions to L^r^JO Although the individual factors cannot be studied in 
isolation, since the coupling between them is due to the finite mass W bosons, one can see 
physical effects of the 'heavier' factors. 

The simplest example is the expectation value of a Wilson loop. At sufficiently long 
distances it is energetically favorable for the heavier string tensions to be screened by W 
bosons. The parameter which controls this is the ratio of string tensions to W masses 
K / m w ~ N 2 m/A, and for any m/A << 1 the almost-A/" = 2 analysis should be valid. 
Thus in the area law regime kL 2 » 1, the Wilson loop will show a crossover from an 
intermediate distance L < mw/ft ~ 1/N 2 m behavior, the sum of terms exp— K n L 2 with 
distinct string tensions, to a long distance behavior governed by the lowest string tension. 

The distinct string tensions are already visible for £77(3). Since the quark charges 
expressed in the monopole basis are (1 0), (—1 1) and (0 — 1), the intermediate range 
fundamental Wilson loop will go as 2 exp —K\L 2 + exp — 2k\L 2 . 

The underlying [/(l)^ -1 symmetry of the effective theory has other physical conse- 
quences as well. It is clearly visible in the light spectrum, the 'glueballs' described by 
the fields of our effective Lagrangian. These are weakly coupled because is small (for 
m « A), and the different £7(1) factors are derivatively coupled (again, controlled by the 
W masses). 

Heavier confined states will vary in each factor as well. Let us consider a qq meson 
state. (We assume for illustrative purposes that the theory with quarks is qualitatively 
similar, which will be true for heavy enough quarks. It is even true for light quarks in 
some cases, for example SU(2) with one flavor. Of course one could talk about states 
containing gluinos in the present theory.) For each single qq state found in a conventional 
analysis (for example the strong coupling expansion), it appears that we would find N 
states, distinguished by the charge of the quark, and bound with different string tensions. 
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The SU(N) gauge theory broken to {/(l)-^ -1 theory has an 5jy discrete gauge sym- 
metry, and one might have thought that this symmetry would somehow remove the multi- 
plicity. However, we found that in the model under study, the dynamics picked a vacuum 
((pi) which spontaneously breaks the symmetry. This affects the expectation values 
and the superpotential (|2.17| ), and thus the string tensions and the physically observable 



'glueball' masses fl2.2Q ) all break the symmetry. The masses of the N meson states contain 
contributions from both the string tensions and the quark masses m q + a^, and (except 
possibly at special values of m, m q and N), there are no degeneracies between them. This 
leaves no doubt that they are distinct physical states. 

The N states are not distinguished by any conserved quantum numbers and thus the 
heavier states are unstable to decay into the lightest. This is mediated by the derivative 



couplings between the U{1) factors and thus a decay amplitude will be A ~ (AE) 2 /m 2 
From the string tensions we can estimate AE ~ V mA and thus the decay rate A ~ m/A 
is controlled by the same arbitrarily small parameter which allowed us to see the distinct 
string tensions. 

The structure of iV-fold split multiplets is a rather unexpected difference between the 
physics of 'almost-TV = 2' supersymmetric gauge theory and non-supersymmetric gauge 
theory. Although definitive results for pure M = 1 gauge theory (or large m in the present 
theory) do not exist, we know of no evidence that this structure would persist, and believe 
it does not. Rather, we believe it is associated with the special features of almost-A/" = 2 
supersymmetry, in particular the light scalar. 

One indication of this is that the natural gauge-invariant operators which create these 
states distinguish the U{1) factors by using the scalar vev. They are the chiral superfields 

O k = QP k (cP)Q. (3.1) 

The polynomials P k {4>) = Cfc 11^(0 - 00 = 2c k T N (</>/2)/(<l> - (p k ) satisfy P k ((f>i) = S ik 
and pick out a single quark charge. 

We might try to identify the split multiplets as bound states with the light scalar. This 
is clearly an appropriate description for scalar mass m >> A. Thus we might hypothesize 
that as we increase m, a split multiplet evolves smoothly into a tower of bound states with 
mass splittings of 0{m) and decay rates 0(1). 

It still may appear that there is a distinction between the picture provided by spon- 
taneously broken [/(l)^ -1 x gauge theory, and the more conventional confining de- 
scription, which allows only [/(l)^ -1 x Sn singlet states. However, there is no invariant 
order parameter distinguishing the two phases, the Higgs phase with Sn spontaneously 
broken, versus the confining phase. This is true even in the theory containing fundamental 
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matter, because the center of Sn is trivial. In other words, the eigenvalues of the adjoint 
scalar themselves transform under the iV (or fundamental) of Sn and every operator in 
the Higgs phase of the theory will correspond to an 5jv singlet operator constructed by 
using the scalar. 

Which description is physically more appropriate depends on the strength of fluctua- 
tions in the Sn sector. In the language of spontaneously broken gauge theory, a configu- 
ration contains domains distinguished by different orderings of the (<fii), and separated by 
domain walls. The discrete gauge symmetry has the consequence that the domain walls 
can end on additional string solutions, with energy scale set by the symmetry breaking 
cii — cij. The language of spontaneously broken gauge theory is appropriate at energies 
low compared to this scale, while at higher energies strings and domain walls can be cre- 
ated freely, fluctuations in the Sn are unsuppressed, and the confining language is more 
appropriate. 

The condition for weak coupling in the Sn sector is the same condition, E « mw, 
that we used earlier and motivated on general grounds of validity of the effective La- 
grangian. 

Baryons will also exist in the theory and it is amusing to note that the [/(l)^ -1 
charges of the quarks are such that the flux tubes would form a chain running from the 
i'th quark to the i + Vst quark in series. This can be contrasted to other pictures proposed 
in the past such as a flux tube with a 'Y' junction. Assuming the W bosons exist as stable 
particles, such a state also exists but is heavier for small m. 

4. The large N limit. 

One of the original motivations for this work was to examine the iV — > oo limit of 
the theory, test the assumptions made in previous work, and evaluate the hypothesis that 
large iV gauge theory is simpler than finite N. One question we need to address is the 
strength of the coupling in the theory. Simple probes of this are connected expectation 
values of the gauge invariant field strength tr F 2 . Standard large N counting predicts 
(tr F 2 tr F 2 ) c < 0(1), (tr F 2 tr F 2 tr F 2 ) c < 0(1/N) and so on. If we calculate 
these expectation values at very long distance we can use the effective Lagrangian ( |2.2| ) to 
evaluate them. The field strength will be determined by its abelian parts. 

Let us first examine the naive semiclassical regime where the electric couplings are 
weak and conventional perturbation theory is valid. To leading order the calculation is 
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governed by the eigenvalues of the kinetic term coefficient r^ . Call these eigenvalues Tj. It 
is not hard to see that in this regime Tj > O(N). The two point function is given by: 



This clearly satisfies large iV counting. But this calculation also illustrates a natural 
mechanism for its breakdown; namely, the possibility of any T{ getting small. But this is 
just what happens near points where monopoles get massless, and so we expect large N 
counting to fail there. Near these points the dual coupling is weak, so perturbation theory 
for these quantities and hence ( |4.1| ) is again reliable. This appears to be an example, of 
which several are already known, of the lack of commutativity of the large iV and infrared 
limits of a theory. It is the hope in ordinary confining gauge theories that these limits will 
in fact commute. 

We now can ask what is the domain of validity of the naive semiclassical regime. For 
concreteness we imagine moving on a line in the moduli space that smoothly connects 
the extreme semiclassical regime with the massless monopole point Cq. (Such paths are 
discussed in more detail below.) We can look for the breakdown of the semiclassical regime 
by examining when the one loop approximation becomes comparable to its first nontrivial 
correction, the one instanton term which is of order A 2N . 

Let us estimate the one instanton correction to aDm = f A. The contour a m 

surrounds the cut connecting the branch points x^m and X2m+i- Semiclassically, i.e., for 

small A, each pair of branch points surrounded by a 7 cycle are located very close to a zero 

of P(x), xzi-i, X2i ~ 4>i- To the required accuracy we can ignore the x dependence in every 

factor of P except (x — <f> m ) and (x — 4>m+i), replacing x in these other factors by either 

4> m or 4>m+i- The difference in results provides a measure of the accuracy of the estimate. 

Choosing we write P(x) ~ \Q(x - </>;)(x - 4>i+i) where Q = Ylj^ i;i+1 {<pi - 4>j)- Using 
A = _L g 9P(x) ^ we can wr ^ e 

2m y ox 



This is closely related to the SU(2) expression for ar>m with the matching condition 
^su(2) = A 2A 7Q 2 . The SU(2) one instanton amplitude is a% m ~ A 4 /(</Wi - </>m) 3 
So a Dm is given by (assuming |0 m+ i - <fi m \ » A) 




(4.1) 




(4.2) 




A 2N 



(4.3) 



Q 2 ((pm+1 
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In the semiclassical region we have taken 0y = v{i — (N + l)/2)8ij. Inserting this in 

I we find roughly 

\ 2N A 

a 1 ' (^_)2^ (44) 

Dm (Nl) 2 v 2N K Nv } 1 } 

As N — > oo, a)/ m — > 0. Note this is not merely the e~ N suppression of instantons. 
Such a suppression alone (reflected in the factor A 2N ) would allow this term to become 
important for v ~ A. The extra AH makes this term negligible until v ~ 1/N. (Of 
course this estimate is too crude to detect possible power of N prefactors of a]^ m .) More 
generally, one instanton corrections are negligible until the splittings between neighboring 
(pi are ~ 1/N. Physically, this suppression arises because of the iV different scales that 
are necessary to give all charged particles at least 0(N ) masses. Instanton contributions 
can be understood as associated with specific SU (2) subgroups, and all of the mass scales 
appear in the matching condition determining the coupling in a given SU(2) subgroup. 
Thus the one loop expression for a Dm remains exact except on a vanishingly small region 
of moduli space as N — > oo. 

Of course this 'vanishingly small' region includes Co and hence is of vital interest 
for our study of confinement. So we now turn to the large iV physics in the immediate 
neighborhood of Co- A striking feature we have already alluded to is the emergence of a 
very large hierarchy in scales near this point. The BPS mass formula gives the Wij masses 
as rriij = \f2\ai — a,j\ with ai = ^(sin ^ — sin ^7 ) — > 2 cos ^ at large N. The smallest 
mass is for Wi^+i which for large N goes as m^i+i — > 27 ^^ sin j^. This is of order ~ 1/N 
for i near the center of the band, and ~ 1/N 2 for i near the edge of the band. 

The mass formula does not guarantee the existence of these W states. There are 
hyperplanes in moduli space on which particular BPS saturated states become unstable 
to decay, and thus the states need not be present in the spectrum on both sides 0,0. 
When the ratio of the central charges Z(q, h) = atq 1 + aoih % associated with two particles 
becomes real, Z(q' , h')/Z(q, h) E R, the particle with charges (q + q' , /i + h') is neutrally 
stable to decay into particles (g, h) and (q' , h'). We should check that these light Ws are 
are present in the spectrum near Cq. In fact they cannot be present at every point near 
Cq. The light Ws are nonlocal with respect to the massless monopoles M m and their mass 
formulae are multivalued around the massless monopole point, as in [Q.* Thus for each 
light W there must be a hyperplane passing through Cq on which it goes unstable. 

This does not exclude the possibility that there exists a path connecting Cq with 
the semiclassical regime on which the Ws remain stable. We have found such a path, 



We thank N. Seiberg for pointing this out to us. 
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generalizing a similar path in || for 577(2). It simply interpolates between the two regimes 
as follows: 



For N > 5 and t « 1 the double zeroes of P 2 — 1 are split slightly and remain on the real 
axis. If we can show that for every t < 1 the zeroes of P 2 — 1 remain distinct, they must 
remain on the real axis, and the assignment of branch cuts and periods can be held fixed. 
We have shown this numerically for N < 21, and the behavior of the zeroes is sufficiently 
simple that we are convinced that this persists for arbitrary N. The accompanying plot 
shows the motion of the zeroes of P 2 — 1 as a function of t for N = 7. 

Thus the a Dm will remain imaginary and the a m will remain real for all t, it will be 
impossible for electrically charged states to decay into magnetically charged states for any 
t, and the presence of the W's in the semiclassical regime implies the existence of light 
electrically charged states near Cq. We have not yet proven the existence of a particular 
charged state. Because all a m are real along the trajectory, all W's are potentially neutrally 
stable. Nevertheless the lightest W's will exist all along the trajectory. This follows from 
the fact that the ordering a m > a n for all m > n is preserved along the trajectory, for 
which the argument is simply that the a m are real, so that if the ordering changes there 
must be a point in moduli space at which a m = a m +i, but this would imply the existence 
of a massless W and a singularity on moduli space not associated with a degeneration of 
the curve, which is impossible. Given this ordering, the particles W mim +i are necessarily 
stable all along the trajectory. 



<t>i(t) = (l-*)&| Co +t( 



N+l 
2 



— i 



(4.5) 



17 



Although the W bosons clearly exist in the TV = 2 theory, the question for the M = 1 
theory is subtle, because we are working so close to a surface on which they go unstable. 
We cannot prove it with our present techniques, and will assume it is true. 

Our effective Lagrangian is only valid for momenta small compared to the lightest W 
mass, p < 1/N 2 . This is another indication of the noncommutativity of the large iV and 
infrared limits in this theory. There should be a signature of this scale in the quantities 
we have already calculated, since the logarithmic running of the dual magnetic coupling 
at the edge of the band due to monopole loops should only start at scales below this mass. 
In fact we see in ( |2.15| ) precisely this effect. 

There are a number of phenomena occurring near Cq. First, as we have seen, instan- 
ton effects can become important roughly when eigenvalue splittings A0 ~ A/N. This 
indicates the initial breakdown of the semiclassical regime. Also, near Cq the monopoles 
become light, at some point becoming comparable in mass to the lightest W they couple 
to. From ( [2. 1U| ) and the inverse of (A.l) we see that to give the monopoles such masses 
O'Dm ~ (A/N) smnm/N, we need to vary Cq as 5fa ~ (A/N) cos ir(i — h)/N. This is a 
variation of fa which is 0(1/ N) compared to its value at Cq. We see that the massless- 
ness of the monopoles is due to delicate cancellations, a mathematical explanation of the 
appearance of multiple scales. In this region one can study the combined dynamics of 
electrically and magnetically charged particles as both their masses are made arbitrarily 
light with their ratio fixed. The first bit of this physics is the varying logarithmic cutoff of 
( p,15|) . A more detailed discussion will follow in the next section. 

We now turn our attention to the confining phase by applying the M = 2 breaking 
perturbation NmTrA 2 discussed earlier. The light Ws dramatically affect the domain of 
m we can control. They cause flux tube breaking and coupling of the different U(l) sectors 
with no suppression at energy scales above their mass, making the theory unmanageable. 
So we are restricted to regime where the confining scale (e.g., a glueball mass) is lighter 
than the W scale. This requires taking m ~ A/N 4 . Naively one might have expected to be 
able to take m ~ A. We see that as N — ► oo the abelian flux tube picture of confinement 
is valid on a vanishingly small part of the m axis. 

In this region there is no reason to think that three Wilson loop expectation values 
are other than order one, and the glueball couplings described by the effective Lagrangian 
are manifestly so. If conventional large N physics holds in the pure M = 1 SYM theory 
obtained as m — > oo a complicated but smooth rearrangement of the theory must take 
place. There certainly will be complicated dynamics occurring for m > A/A 4 . Whether 
conventional expectations are realized we cannot say from this analysis. It is conceivable 
that order one couplings could persist as m — > oo. 
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If conventional expectations are valid, they would probably hold until string tensions 
become of order light W masses. Since these masses go to zero as N — > oo the region of 
validity of conventional expectations would occupy the whole m axis at N = oo. There 
would be a sharp transition between the coulomb Af = 2 theory described by one loop 
physics and the confining theory described by large N lore. The light abelian monopole 
region that smoothly interpolates between them would have disappeared. 



5. A scaling trajectory. 

So far, our discussion of the large N limit has been largely qualitative. We saw that 
as we move in the Af = 2 moduli space towards Co, many things happen. We start with a 
simple U (l)^ -1 gauge theory in which many calculations are exact at one loop. Gradually 
instantons turn on, and at the same time monopoles become light, eventually crossing 
the W masses to become the lightest degrees of freedom. They will then drive the beta 
function and cause the U(l) gauge factors to decouple in the monopole basis, different 
from the original basis in which ((f)) was diagonal. Our analysis suggested that all of this 
happened in the scaling regime \(f>i — ~ 1/N . 

We will now study the relationship between these phenomena more systematically by 
introducing the following scaling trajectory for the eigenvalues of ((f)): 

Us) = e s/N <t>i\c =e s / N 2cos^H . (5.1) 

Introducing the scaled coordinate z = e~ s / N x we find the polynomial is of the form 
P(x, s) = e s P(z) and the curve is given by y 2 = e 2s (P 2 (z) - e~ 2s ) where P(z) = P(z)\c . 
As for our previous trajectory, the branch point evolution is smooth: taking z = 2 cos 9 we 
have P(z) = cos N9 and the branch points are at cos N9 = ±e~ s . 

The trajectory is entirely in the scaling regime and even with s large the outlying 
\<f>i — <f>i+i\ remain 0(1/N 2 ). However its essential property is that for s large, the length 
of the 7 cycles shrinks to zero (in terms of x, as 2e~ s ) and thus we recover semiclassical 
results. 

The form A becomes 

_ e s / N zdz dP 

1 e s / N cos 9 d(cosN9) 
ni Vcos 2 N9 - e~ 2s 

19 




e 



The large N simplification occurs because the integrals over the a and 7 cycles only cover a 
range of ~ 1 /N. Thus we can expand the cos factor to isolate the leading N dependence. 
For ai and ao m we obtain* (using u = cos NO) 



a, 



e s ' N 2 cos Bi 



TV 



e s ' N 2cos6> 7 - 



du 



-2s 



U 



(5.3) 



and 



7T 



AW 



M=l 



sin 6 r 



du arccos u 



\/u 2 — e~ 



2s 



(5.4) 



We see that we can make monopoles much heavier than W's along this trajectory. 
Monopole masses are comparable to the lightest W they couple to for s ~ 1. One loop 
contributions going like log (a/ A) will show s dependence ~ s/N, and we see from (|5.4j) 
that at leading order in N the monopole and W masses are given exactly by the one loop 
result, all the way to the massless monopole point! (We will see below that this is a special 
property of this trajectory and more generally trajectories 5(f>i ~ cos7r/c(z — \)/N with 
k ~ O(N ).) 



* The integral h(s) is most easily evaluated by using partial integration to establish Ji(s) = 
—Ii(—s). The only solution to this compatible with the analytic properties of I\ is I\ = s. We 
are grateful to Alyosha and Sasha Zamolodchikov for providing this argument. 



20 



To compute r, we might use 



dX 

d(j>i 



,s/N 



sin# 



cos N9 



2n cos 9 - cos d t Vcos 2 N9 - e~ 2s 



d9. 



(5.5) 



There is a subtlety in the large iV expansion of this form, due to the pole in the prefactor 
at 9 = 9i. Although the pole itself is cancelled by a zero of cos N9, it is not justified to 
treat the prefactor as slowly varying. 

We will deal with this by using the basis of variations 



N 



dX dX wk(i 

= 1^ 7HT cos " 



dt k 



N 



(5.6) 



We will find that using this basis for 'electric' quantities and the basis sin k9 m for 'magnetic' 
quantities, the period matrices become simple and the kinetic term diagonal, for all s and 
at finite N. 

Before doing this in detail, let us suggest why this gives simple results. From (A.l) 
and related formulas, it is clearly an appropriate basis for the massless monopole point 
s = 0. On the other hand, the large s limit is expected to reproduce the semiclassical 
results ddi/dtfij = 5ij and daoi/dcpj = d 2 J r /daidaj. From (|2.3|) and ( |5.3[ ) this is 



da 



Di 



T 



(1) 



one loop 



LI 



2tt 
2s 



&ij ^^tik — (1 — 5ij)t 



(5.7) 



= ^ + log (2 cos6>i — 2 cos6>j) 2 . 



The point is that the one-loop kinetic term is also diagonal in the cosine basis, by results 
similar to those in the appendix, making it possible for r along the entire trajectory to be 
diagonal. Let us estimate its eigenvalues in the large limit: the difference — t 

B m =ij +0(1/N) and so (A.l) implies that Uj will have eigenvalues t h 
are corrections to this estimate, of 0(k/N)). 
We proceed to compute the periods of 



2sin7rfc/7V' 



(There 



dt k 



N 



EdX nk(i - i) 
^cos 



d(pi 



N 



i=i 

1 N sm(N — k)9 
2k Vcos 2 N9 - e~ 2s 
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(5.8) 



We have (a = arccos e s , 6 = arcsine s , K = k/N) 



<9t fc 7T 7_ & \fe~ 2s - sin 2 6> 

1 /" & cos /c6*j cos(l — k)6> + sin k9j sin(l — k)9 

n J-b \fe~ 2s - sin 2 9 

= F(k, s) cos k9j 

1 /" & cos(l - k)0 



1 

F(«, s) = - d9 

7T 7-6 



(5.9) 



^ 7-6 Ve -25 - sin 2 
— sin — (— log s) + 1 — O(k) s^O 

TV 2 

1 + -^— : — -e 2s + ... s -> oo 



and similarly 



Bmk(s) = da ® m = G(k, s) sin k§ r 



i r 

G(k, s) = - d6 

7T 7-a 



1 /■» cos(1-k)0' 



Vcos 2 0' — e 2s 



(5.10) 



71 K 

— > 2s sin — + 1 — O(k) + . . . s — > oo. 

We can now compute TD mn (s). Using (A. 5) we see that T£> mn (s) is diagonal in the 
sinkOm basis: 

TDmn{s) = ^ A mk B fc J = ^ q^AjkB^, 

y^TD mn (g) sin/c(9 n = t Dk (s) sin/c/9 r 

n 

with eigenvalues 



(5.11) 



W.) = ^-151^ • (5-12) 
We verify that td k (s) — > — ^ logs, independent of /c, as s — > 0. Similarly 

T ij(s) = J2 h ^ BmkA k J 1 

m,k 

yV^(s) cos k6j = r K (s) coskOi (5.13) 
1 G(k,s) 



2sin^ s) 
22 



These formulae give an exact description of the physics along the scaling trajectory from 
the semi classical domain all the way to Co, even at finite N. 

Note the crucial fact that the eigenvalues t k (s) are not simply the inverses of the 
eigenvalues td k (s)1 The change of basis from electric (ij) to magnetic (mn) couplings 
introduces an extra 4 sin a wk ■ This will have important consequences and seems quite 
mysterious; it may be relevant that (for example) in the magnetic basis, the W bosons 
responsible for the one loop contribution to r have charges with alternating signs (such as 
(1 — 1 —11)), and in summing their loop effects there will be large cancellations. 

Let us consider the leading correction to the logarithm from the monopole beta func- 
tion at small s, which come from the constant terms in F and G: 

td k (s) = ~ log s + + Oil) + .... (5.14) 

2% 2 sin *y 

For k ~ O(N ) (k ~ 1/N) and s ~ O(N ), this is a large correction, of O(N). It 
is a threshold effect and eventually the monopole-induced running of the coupling will 
overwhelm it, but we need to go to scales s ~ e~ N for this to happen. Thus these effects 
are quite important at large N. 

The same constants are responsible for the s-independent term at large s in 

r,(s) = 2s + -^—,+0(l) + .... (5.15) 

— ■ bill 2 

We see by comparison with (|5.7| ) that this is just r one i oop . Rather magically, the low k 
'modes' manage to be simultaneously weak in electric and magnetic variables. 

Since the constant term in td k has the same form, the form of To in the monopole 
basis can be inferred by analogy to r: 

TDmn —5 mn logs+ — log(2 sin 9 m - 2 sin 6 n ) 2 + 0(1) + .... (5.16) 

This matrix is not diagonal and so a 'magnetic observer' would see peculiar transitions 
between the m quanta. This observer might posit the existence of light solitons carrying 
dual charge in four m factors (as mentioned above) that would mediate these transitions. 

The similarity of (|5.16|) to the one-loop electric r is very intriguing and leads us to 
wonder if there is any sense in which it is a one-loop effect. The theory is certainly much 
more 'self-dual' than one would have guessed. 

We now plot the scaling functions to illustrate their structure. Note that the — ^- log s 
behavior characteristic of the monopole one loop beta function sets in around s = 1 
uniformly in k. This is just at the energy scale set by the inverse 'size' of the monopole, as 
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discussed earlier. The large residual s independent term in ( |5.14j ) is responsible for the k 
dependent offsets of the curves. It is clear from the graph that one must go to s ~ e~ l / K 
before this offset is negligible. 

The off diagonal terms in ( 5.16 ) influence the energies of states in the confining phase 
obtained by adding the M = 2 breaking perturbation NmTrA 2 . The basic structure 
of the flux tubes is determined by the monopole condensate which is determined from 
the superpotential alone and hence does not depend on 7"D mn . It is still given by ( |2.19|) . 
The string tensions do depend on Tomn and will be affected by the mixings described by 
( |5.16|) . We have not yet studied this problem, but we have examined the closely related 
problem of the dual photon masses, by numerically diagonalizing the mass matrix. We 
found that the quantitative relation ( |2.2U[ ) is altered but that the scales of the masses 
remain unchanged. We expect the same for the string tensions. Note that these changes 
are significant for m ~ e~ N , far below the energies where string breaking due to light W 
pairs occur (m ~ 1/N ). 

We remind the reader that some quantities (a, and a Dm) were exact at one loop along 
the scaling trajectory. This occurred because the corrections to the one-loop results for 
their derivatives, from (|5.9| )-( |5TT0| ), were suppressed by powers of k, so would happen for 
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any trajectory with k ~ 1/N. If the same thing happened in other large A" (supersym- 
metric) theories, for even a single trajectory into the true ground state (here the massless 
monopole point) , it would have the important consequence that the ground state could be 
found just knowing the one-loop effective action. 



6. Conclusions 

Using the solution of Af = 2 supersymmetric SU (N) gauge theory of [0,|3],|[] , we have 
given an overview of some of the physics visible in the low energy effective Lagrangian. In 
particular we studied the Af = 1 theory obtained by giving a mass to the chiral superfield, 
and extended Seiberg and Witten's explicit derivation of a monopole condensation model 
of confinement for SU(2) to SU(N). 

The Lagrangian relevant for this confining theory is a weakly coupled {/(l)^ -1 gauge 
theory with one light monopole hypermultiplet in each factor. The Sn discrete gauge sym- 
metry that permutes the U(l) factors is spontaneously broken. This phenomenon pervades 
the physics of this system. In particular, there is a spectrum of different string tensions in 
the different factors, K n ~ mAN sin ^2 anc i a spectrum of different light W boson masses 
~ -4 sin that connect the factors. This differs in many ways from our expectations 
about the pure Af = 1 theory. It is possible for a dramatic but smooth rearrangement to 
take place for m ~ A that will connect these two descriptions. Perhaps this rearrangement 
can be understood as a smooth crossover from a regime of spontaneously broken Sn gauge 
symmetry to one where the Sn gauge symmetry is 'confined.' 

We studied the large iV limit in detail. The hierarchy of scales present for finite N 
becomes very large, as the lightest W mass ~ A/N 2 . We found one signature of these very 
light particles in the monopole 'size' which controls the energy scale of onset of pertur- 
bative monopole coupling constant renormalization. Another very surprising signature is 
the off-diagonal terms in (|5.16|) which cause a coupling of the different magnetic factors. 



Understanding this phenomenon from the magnetic point of view, presumably as the result 
of light electric solitons in a weakly coupled magnetic theory, would be very interesting. 

At large A" the one loop result becomes exact almost everywhere on the Af = 2 moduli 
space. We identified a scaling regime very close (~ 1/N) to the massless monopole point 
Co where instanton and monopole effects survive the large A" limit and were able to give 
simple exact formulas for monopole and W masses and coupling constants along a scaling 
trajectory through this regime. 

The coupling almost everywhere on the Af = 2 moduli space is 1/N, as expected, but 
near Cq becomes large, contradicting standard large A" lore. The infrared divergence of the 
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effective U(l) couplings, which does not commute with the large iV limit, is the source of 
this. Such a phenomenon would seem to apply to any continuous monopole condensation 
picture of confinement, supersymmetric or otherwise. The region of this violation goes to 
zero as N — > oo. 

In the confining 'almost Af = 2' theory, the abelian flux tube picture breaks down 
because of W pair creation when confining scales become comparable to W scales. This 
occurs for m ~ A/A 4 , i.e., a vanishingly small region of the m axis. Conventional large N 
lore could be recovered in a smooth, complicated crossover we cannot control. As N — > oo 
the region described by nontrivial light monopole physics shrinks to zero, possibly leaving 
an abrupt transition between one loop physics and large N pure N = 1 SYM behavior. 

Some features of this model are reminiscent of standard hadronic phenomenology 
while others look very different (such as the weakly coupled glueballs), but overall we find 
the results encouraging for the idea of using supersymmetric gauge theory as a solvable 
starting point for attempts to model QCD dynamics. Surely interesting supersymmetric 
analogs of many other problems of strongly coupled gauge theory and QCD physics can 
be found, and we believe it will be very fruitful to study them using these techniques. 
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Appendix A. Trigonometric identities. 

In section 2 we use 



N-l 




(A.l) 



m=l 



which can be checked by writing (z. 



rn 




Wi = e l6i ) 




(A.2) 
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We convert the sum to run from to 2N — 1, for which e n ' ikrn ' N = 2N5 k ^ Q ( rnod2 N)i by 
taking a sum ^m=i involving the second term in ( |A.2| ) and taking m — > 2N — m. Then 

N-l , ^ 2N-1 _i 

71KTYI 1 \ W- 



-iirkm / N p iirkm/ N 



N 2N ^ wr'-zr, '' 



w 2 t N ~ k (A.3) 



1 - w1 N 1 - wf N 
2 ' 

vZV-1 



Using ( |A.l] ), the orthogonality of the basis smnkm/N under Yu m =i' an< ^ the orthog- 
onality of the basis costt/c(z — \)/N under 52 i=1 , we also have 



B mi cos /c6*i = sin k9 m . (A. 4) 



In section 5 we use the change of basis 

<7™ s i n k9 m = sin k9 m= i — sin k9 m= i-\ = 2 sin cos (A. 5) 

771 

We next evaluate the sum 

H9) = if Sind cos fc^ = Sln(iV ~ fe)g . (A.6) 

First, for fe = 0, the sum is 

^^ lo sIl^ os ^ ~ cos #) = ^^logP(x) =tanm (A.7) 

3 

For k integer, the sum is a periodic function of 9 with residues cos k9 at the points 9 — 9j, 
so it must take the form 

f k (9)= tan N9 cos fc0 + finite. ( A.8) 

From the above, the finite part is known at the points 9 = 9 m to be — sin k9. This is true 
of the function 

f k (9)= tan N9 cos k9 - sin k9 ( A.9) 

(equal to ( |A.(j| )) plus any periodic function vanishing at all of the 9 m , such as sin2A^6>. 
Such terms can be excluded by considering the derivative d/d9 at 9 = 0. 
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